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ABSTRACT 

Context. GJ436b is the first extrasolar planet discovered that resembles Neptune in mass and radius. Two more are known (HAT-P- lib 
and Kepler-4b), and many more are expected to be found in the upcoming years. The particularly interesting property of Neptune- 
sized planets is that their mass M p and radius R p are close to theoretical M-R relations of water planets. Given M p , R p , and equilibrium 
temperature, however, various internal compositions are possible. 

Aims. A broad set of interior structure models is presented here that illustrates the dependence of internal composition and possible 
phases of water occurring in presumably water-rich planets, such as GJ 436b on the uncertainty in atmospheric temperature profile 
and mean density. We show how the set of solutions can be narrowed down if theoretical constraints from formation and model 
atmospheres are applied or potentially observational constraints for the atmospheric metallicity Z\ and the tidal Love number k 2 . 
Methods. We model the interior by assuming either three layers (hydrogen-helium envelope, water layer, rock core) or two layers 
(H/He/ILO envelope, rocky core). For water, we use the equation of state ILO-REOS based on finite temperature - density functional 
theory - molecular dynamics (FT-DFT-MD) simulations. 

Results. Some admixture of H/He appears mandatory for explaining the measured radius. For the warmest considered models, the 
H/He mass fraction can reduce to 10~ 3 , still extending over ~ 0.7R e . If water occurs, it will be essentially in the plasma phase or in the 
superionic phase, but not in an ice phase. Metal-free envelope models have 0.02 <k 2 < 0.2, and the core mass cannot be determined 
from a measurement of k 2 . In contrast, models with 0.3 < ki < 0.82 require high metallicities Z\ < 0.89 in the outer envelope. The 
uncertainty in core mass decreases to 0.4M p , if k 2 > 0.3, and further to 0.2M p , if k 2 > 0.5, and core mass and Z\ become sensitive 
functions of k 2 . 

Conclusions. To further narrow the set of solutions, a proper treatment of the atmosphere and the evolution is necessary. We encourage 
efforts to observationally determine the atmospheric metallicity and the Love number k 2 . 

Key words, planets and satellites: interiors - planets and satellites: individual: GJ 436b - EOS: water equation of state 



1. Introduction 



The plan et GJ436b dButler et al.l2004 has attracted much atten- 
tion (e .g. iManess et al. I l2007t iBaraffe et alj|2008t IBatygin et al. 
I2009al) because it is the first Nept une-mass extrasolar planet ob- 
served transiting its parent star dGillon et alJ l2007bf) . Follow- 
up Doppler observations and infrared p hotometry dat a from 
primary transit and secondary eclipse (iDemorv et al.l l2007t 
iDeming et al.i r2007) have allowed its mass M p , radius R p , and 
day-side brightness temperature to be determined to within 14%, 
13% and 25%, respectively. For Neptune-size extrasolar planets, 
these parameters permit a variety of possible compositions rang- 
ing from a rock-iron planet surrounded by a gaseo us H/He layer 
to a planet that is substantia lly composed of water (lAdams et alJ 
l2008t iFortnev et aT]l2007al) . For instance, for the Neptune-size 
solar planets Uranus and Neptune, more observational con- 
straints such as the gravitational moments are available, but their 
bulk composition is still unknown due to the a mbiguity of the 
mean density of water and H/He/rock mixtures dPodolak et alJ 

I2000h. 

iRogers & Seagert d2010ah quantify the uncertainty in the 
composition of Neptune-mass planets by taking the uncertainties 
in M p , R p , and intrinsic luminosity into account. For GJ 436b, 
they obtain mass fractions of H/He, water, and rocks, where 
rocks comprise silicates and iron, of 3.6-14.5%, 0-96.4%, 
and 4-90%, respectively. They find a particularly strong de- 



pendence of compos ition on the atmospheric thermal profile. 
iFigueira et al.l d2009l) are able to further constrain the compo- 
sition by comparing with possible formation histories. Their 
formation models exclude water-less models and predict initial 
mass fractions of H/He, ice, and rock of 10-20%, 45-70%, and 
17-40%, respectively. This demonstrates the necessity of fur- 
ther constraints to distinguish between various interior solutions. 
Su ch an additional observable would be the tidal Love number 
k 2 dRagozzine & Wolfll2009H"Batvgin et al.ll2009bl). For the two- 
planet system HAT-PI 3b,c . iBatygin et alJ d2009bh show how 
high-precision measurements of the orbital and planetary param- 
eters can translate into information about the planetary interior 
of the inner planet, such as parametrized in terms of a core mass. 
GJ 436b's high ecce ntricity suggests th e presence of such an ad- 
ditional companion dRibas et al.ll2008l). Predictions for the per- 
turbers's mass of « 10M m dCoughlin et al]|2008l: IBatygin et al.l 
l2009al) and eccent ricity are just within the scope of current de- 
tection thresholds dBatvgin et al.ll2009ah . so the value of the ob- 
servable kn could be inferred. In this paper we discuss the infor- 
mation content of k 2 with respect to the core mass and metallicity 
of GJ 436b, and its implication for a secondary planet. 

Core mass and metallicity are calculated for a broad range of 
models that aim to embrace the full set of solutions that are al- 
lowed by the uncertainty in M p , R p and the atmospheric temper- 
ature profile. The material components that GJ 436b is assumed 



2 



Nettelmann et al.: Interior models of GJ 436b 



to be made of are rocks, confined to a rocky core, water, whether 
confined to an inner envelope or uniformly mixed into an outer 
hydrogen-helium envelope, and hydrogen and helium. For H, 
He, and water, we apply the Linear Mixing Rostock Equatio n 
Of State (LM-REOS) described in INettelmann et al.l (120081) . 
which is based on finite temperature - density functional theory 
- molecular dynamics (FT-DFT-M D) simulations for the compo- 
nents H, He, and H2O, (s ee iHolst et ail 120081: iKietzmann et alj 
l2007t iFrench" et al.ll2009l) . In particular, the phase diagram of 
water has been calculated recen tly up to pressures of 100 Mbar 
and temperatures of 40000 K (IFrench et alj|2009l) . so that we 
can derive the possible phases of water in presumably water-rich 
planets such as GJ436b in dependence on the uncertainties men- 
tioned above. 

We describe and discuss the observational constraints ap- 
plied in § 12.11 the equation of state data in § 12.21 and our method 
of generating structure models in § 12.31 The range of resulting 
compositions are presented and discussed in § 13.11 in the con- 
text of the phase diagram of water and constraints from forma- 
tion theory and model atmospheres. The range of resulting core 
masses and metallicities is presented and discussed in § 13.21 to- 
gether with the potentially observable Love number Ic2- We find 
that, at deep interior temperatures and pressures in GJ 436b, 
finite-temperature effects of the water EOS are important. In 
§ !3.3l we present new M-R relations of warm water planets. Our 
conclusions are discussed and summarized in §3] 

2. Input data and model construction 

2.1. Observational constraints 

Observational data of the mass M p , radius R p , and surface tem- 
perature of GJ 436b have been provided by serveral authors as 
summarized below. It is interesting to note that different deter- 
minations of M p and R p arise mostly from different assumptions 
ab out the stellar mass M , and radius R*. 

iDeming et alJ d2007l) determined a planet-to-star radius ratio 
R P IR, of 0.0839 + 0.0005 from the depth in the precise Spitzer 
transit photometry light curve. They then used the Spitzer light 
curve and radial velocity data to find the best fit to the stellar 
radius versus the stellar mass and compared this relation to the 
empirical M-R relation M» = R t to find M, = R t = 0.47 ± 0.02 
(in solar units) and consequently a planet radius R p = 4.33+0.18 
Earth radii (R®). From the amplitude of the secondary eclipse 
observed with Spitzer, they also determined a day-side bright- 
ness temperature of the planet of 712 ± 36 K. They also give 
M p = 0.070+0.003 Jupiter masses (M s ) as best fit t o all observa- 
tions. A ssumi ng M« = 0.44M H and R* = 0.463^. iGillon et all 
d2007albt) and lDemorv et al.l d2007l) find similar planetary values 
wit hin the error bars by analyzing the Spitzer transit data. 

iTorres et"al] (120081) use more recent results from stellar evo- 
lution models that match the observed stellar luminosity derived 
from spectroscopic measurements of the effective stellar temper- 
ature and the mean density derived from photometric measure- 
ments of a//?*, where a is the planet-star separation. They obtain 
as best values M, = 0.452 + 0.014M o and /?» = 0.464 ± 0.01 1R 
for GJ 436 with error bars covering previous estimates, and 
M p = 0.0729 + 0.0025Mj, andR p = 0.37 67 + °™fo, for GJ4 36b, 
whe re Rj is the equatoria l Jupiter radius. iTorres et all d2008h and 
also lDeming et al.l (12007) assume the effective temperature r, >e ff 
of GJ 436 to be 3350 + 300 K, which we think is too low and not 
constrained enough for an M3V star with a V-K color index of 
4.7 mag (V=10.68, Simbad; K=6.073, 2M ASS). For such an M 
star, we obtain T. t e g = 3470 + 100 K from iKenvon & Hartmanl 



(119951) in agreement with iBean et all d2006l) . who compared 
their optical high-resolution spectra of GJ 436 with the lat- 
est PHOENIX model atmospheres to obtain 7\ eff = 3480 K. 
This leads to a slightly higher planet equilibrium temperature 
Teg = T* ieS (RJ2a) l/2 of 673 ± 20 K instead of 649 ± 60 K, 
assuming zero-albedo. Within the error bars, however, the equi- 
librium temperatures agree with each other and also with the 
plane t's effective t emperature (7y e ff = 712±36K. lDeming et al.l 
120071: 717 + 35 K. lDemorv et al.ll2007l) obtained from Spitzer in- 
frared data, but no longer does so if assuming a non-zero albedo 
a p = 0.3 as of Uranus yielding T eq = 535 K, which would then 
imply a significant intrinsic luminosity. 



IBean et ail (120081) have repeatedly observed the transit with 
the Hubble Space Telescope (HST) Fine Guidance Sensor (FGS) 
and fitted the transit light curves to directly obtain the stellar ra- 
dius, as well as the planetary radius, orbital period, and inclina- 
tion. They find R p = 4.90+q 'iiR®, i.e. larger than all other radii 
measured by others. Our baseline m odels have M„ - 2 3.17M® 
and R p = 4.22/? ffi in accordance with lTorres et al ] (120081) . 



2.2. EOS 

Interior models presented in this work are composed of H, 
He, H2O, and rocks. The rock -EOS used here is an an alyti- 
cal pressure-density relation by iHubbard & Marlevl d 1989h that 
approximates a theoretical EOS of a mixture of 38% Si02, 
25% MgO, 25% FeS, and 12% FeO appropriate for Jovian core 
conditions, i.e. for T ~ 10 4 K. This rock-EOS has been applied 
extensively to the core region of the solar sy stem giant planets 



'eg 

dHubbard & Marlevl 1989tlGuillot et al.ll 19941) 

Mixtures of H, He, and H2O are generated by linear mix- 
ing. For these components we use H-, He-, and H2O-REOS, re- 
spectively, which are based on FT-DFT-MD simulations at those 
densities where correlation effects are import ant and on chem- 
ical m odels in complementary regions, see INettelmann et al.l 
(120081) and references therein for a more detailed description. 
In particular, H2O-REOS is a combination of different finite- 
temperature EOS covering phases ice I and liquid water by 
fits to accu rate experimental data , water vapor using Sesame 
EOS 7150 dLvon & Johnsoniri992l) . supercritical molecular wa- 
ter (Sesame 7150 and FT-DFT-MD), fractions of ice VII and X, 
and large regions of superionic water and water plasma up to 
20 g cm 3 and 40000 K (FT-DFT-MD). All interior models pre- 
sented here are too warm for ice phases to occur. While real 
planets may have accreted water in combination with methane 
and ammonia (together conveniently named ices regardless of 
their thermodynamic phase at planetary interior conditions), we 
consider the ab initio EOS H2O-REOS an improvement over for- 
mer ;ce-EOS, because temperature effects are properly taken into 
account. 



For comparison, lAdams et al.l ( 2008b and iRogers & Seageri 
d2010al) use EOS data bv lSeager et alJ (120071) for their GJ 436b 
models, where the iron and silicate EOS reproduce experimen- 
tal data at P < 2 Mbar and the water EOS is in a ddition based 
on DF T-MD data for water-ice phases VIII and X. iBaraffe et al.l 
d2008l) extensively studied the effect of different available water 
EOS, i.e. of Sesame and ANEOS water EOS. Isotherms of these 
EOS show a significantly more compres sible behavior in th e 
Mbar region than the new H2O-REOS, see IFrench et al.l d2009l) . 
thereby requiring a higher mass fraction of the light H/He com- 
ponent to yield the same planetary mean density. 
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2.3. Structure assumptions 

We obtained interior models by integrating the equation of 
hydrostatic equilibrium of a gravitating sphere dP/dr = 
-Gm(r)p(r) / r 2 and the equation of mass conservation dm/dr = 
4nr 2 p(r) inward, where P is pressure, p mass density, m the 
mass, and r the radial coordinate. These equations require three 
conditions at the outer boundary r\ = R p for the variables m, P, 
and p. Clearly, m(r\) = M p . In accordan ce with prior int erior 
structure models of solar giant planets (e.g. lStevensonll982T) . we 
define the surface as the 1-bar pressure level, thus P{r\) = 1 bar. 
For pi = p{P\, T\) from the EOS we need the thermodynamical 
surface temperature T\ which we consider a variable parame- 
ter. Since a real planet has to obey the inner boundary condition 
m(0) = 0, we make the additional assumption that a core ex- 
ists and choose the respective core mass M core such that total 
mass conservation is ensured. We define a layer by its homoge- 
neous composition, where pressure P„-(„+i) and temperature at 
the transition from layer No. n to No. n + 1 are continuous, while 
density and entropy change discontinuously. 

In this work we consider two-layer (2L) and three-layer (3L) 
models. For our 3L models, the innermost layer No. 3 is a core 
of rocks, layer No. 2 is an adiabatic inner envelope of water 
with envelope metallicity Z% = 1, and layer No. 1 is an outer 
envelope of linearly mixed hydrogen and helium with an He 
mass fraction Y = Mu e /(Mn + Mh c ) = 0.27 as is character- 
istic f or the protosolar clo ud and envelope metallicity Z\ = 0. 
As in lAdams et al. I (l2008h . we assume an adiabatic profile for 
P > 1 kbar that is located in the outer envelope. This struc- 
tur e type is in line with pioneering Uranus and Neptune models 
by iHubbard & MacFarland d 1980b . who investigated a layered 
interior differentiated into a central rocky (silicates+iron) core 
overlaid by an ice shell (H2O, CH4, NH3) and surrounded by a 
H/He envelope. While such models do not reprod uce the gravita- 
tional moments J 2 and J 4 of Uranus and Neptune dHubbard et al.l 
Il995l) . implying an admixture of heavy elements into the outer 
H/He envelope, as well as some admixture of H/He into the in- 
ner, ice-rich envelope, they are a reasonable starting point when 
additional gravity field data are not at hand. Given only M p , R p , 
and 7^ e ff and a location of GJ436b close to theoretical M-R rela- 
tions for pure water planets, we expect a broad range of solutions 
even if observational uncertainties in M p and R p are neglegted. 
For our 3L structure type models, we vary the transition pressure 
Pi _2 from the outer to the inner envelope between the bounds as 
defined by zero-mass water layer and zero-mass rocky core. 

Alternatively, we calculated two-layer models where water is 
homogeneously mixed into the H/He envelope. For this structure 
type, we varied the envelope metallicity Z\ := MH,o.env/^env be- 
tween the bounds as defined by Z\ =0 and a zero-mass core on 
the other hand. These two kinds of simple structure assumptions 
bracket the uncertainty regarding the distribution of ices within 
the envelopes. The presence of heavier elements such as silicates 
in the H/He/H20 envelope or in the water layer of the 3L models 
would act to reduce the amount of water in favor of the amount 
of H/He needed to reproduce the radius. These models therefore 
give a lower limit for the H/He mass fraction in GJ 436b. 

An important question is the temperature profile between 1 
and 1000 bar, where atmospheric models for irradiated planets 
based on non-gray radiative-convective equ ilibrium calculations 
dFortnev et alj|2007at iBurrows etal] 120081) predict an isother- 
mal region. Adiabatic atmospheres, on the other hand, are gen- 
erally assumed for all solar giant planets owing to weak irra- 
diation and strong molecular absorption in that pressure region 
dGuillot et al.ll 19941) making energy transport by radiation ineffi- 



cient. For Jupiter, this assumption is confirmed by Galileo probe 
entry measurements and for Uranus & Ne ptune by Voyager in - 
frared and radio occultation observations dGautier et al.ll 1995b . 
Since we aim to not narrow the possible set of solutions in ad- 
vance by the choice of a specific model atmosphere, we cover 
these limits in the case of our 3L models by both calculat- 
ing isothermal and adiabatic models in the pressure range 1 to 
1000 bar representing the extreme cases of strong or no influ- 
ence of irradiation. For both 2L and 3L models, we vary T\ 
between 300 K and 1300 K. Cold interior models of GJ 436b, 
in particular with water cores at 300 K overlaid by an adia- 
ba tic, hence an even co lder, H/He atmosphere have been studied 
by Baraffe et al. d2008l) with respect to the evolution timescale. 
ISpiegel et al. ( 20101) . on the other hand, find isothermal atmo- 
spheres at 1300 K and 1-100 bar, independent of the day-to 
nightside energy redistribution in the atmosphere of GJ436b. For 
our 2L models, we also vary T\ := 7(1 bar) from 300 to 1300 K 
but only consider temperature profiles with T\ = 7ioo constant, 
where 7ioo := 7X100 bar), and adiabatic for P > 100 bar. 

We subdivide our structure type models into 1 1 series labeled 
A, B,. . . , K representing the uncertainties in T\, 7iooo, and mean 
density p p , where 7iooo := 7X1000 bar). For series A-I, we take 
M p = 0.0729Mj = 23.17M ffi and R p = 0.3767flj = 4.22fl e ; for 
series J, we take M p = 21.3M® and R p = 5.35P e corresponding 
to the minimum value of p p ; for series K, we take M p = 24M e 
and R p = 47?$ representing a maximal p p . They are modifica- 
tions of series E intended to investigate the effect of the mean 
density. The input characteristics of our model series A-K are 
listed in Table Q] 



Table 1. Input characteristics of model series A-K. 



series 


r, 


v r 


M p 


R P 


N 




[K] 




W B ] 




(layers) 


A 


300 


adiabatic 


23.17 


4.22 


3 


B 


300 


isothermal 1-10 3 bar 


23.17 


4.22 


3 


C 


300 


isothermal 1-10 2 bar 


23.17 


4.22 


2 


D 


700 


adiabatic 


23.17 


4.22 


3 


B 


700 


isothermal 1-10 3 bar 


23.17 


4.22 


3 


F 


700 


isothermal 1-10 2 bar 


23.17 


4.22 


2 


G 


1300 


adiabatic 


23.17 


4.22 


3 


H 


1300 


isothermal 1-10 3 bar 


23.17 


4.22 


3 


1 


1300 


isothermal 1-10 2 bar 


23.17 


4.22 


2 


J 


700 


isothermal 1-10 3 bar 


21.3 


5.35 


3 


K 


700 


isothermal 1-10 3 bar 


24.0 


4.0 


3 


"or series A-I, M„ 


and R„ are taken from Torres et 


al. (2008), yieldin 



a mean density p p = 1.69 g cm , while in series J, p p = 0.45 g cm 
and in series K, p p = 2.22 g cirT 3 . 



2.3.1 . Love number k 2 

For selected models, we calculated the Love number k 2 . 
This planetary property quantifies the strength of the planet's 
quadrupolic gravity field deformation in response to an ex- 
ternal massive body. In case of the system GJ 436 and 
GJ 436b, M t causes a tide-raising potential W(r) — 2w=2 ^« — 
{GM t /a)Y Jn= 2( r / a )"Pn(cosi//), where a is the (time-dependent) 
distance of the center of masses, r the radial coordinate of a 
point inside the planet in a planet-centered coordinate system, 
ip the angle between a planetary mass element at r and M, at 
a, and P„ are Legendre polynomials. Each external potential's 
pole moment W„ induces a disturbance V„ = K n (r)W„ of the 
planetary potential, with k„ = K n (R p ) taken at the surface R p the 
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Love numb ers. Relevant equations to c alculate Love numbers 
are given in lZharkov & Trubitsvnld 19781 §48). 

As the gravitational moments Jj^, the Love numbers are 
uniquely defined by the planet's internal density distribution. To 
first order in the exp ansion of a plan et's potential, Ic2 is propor- 
tional to Jx (see e.g. iHubbardll 19841) . so that measuring ki pro- 
vides us with an additional constraint equivalent to J2 for the 
solar system giant planets. In particular, £2 is -as is J 2- most 
sensitive to the central condensation of a planet, as parametrized 
adequately by the mass of the core within a two-layer approach. 
However, as J2 and £2 depend on the accumulated density dis- 
tribution p(r), the inverse process of deriving p(r) from J2 or 
£2 is not unique. Even in the most fortunate case of enough 
moments being available to accurately determine the internal 
density distribution -which i s still uncertain even for Jupiter 
dFortnev & Nettelmann! l2010h - we have to face the ambiguity 
from various materials yielding the same density when mixed. 
For example, the common perception of Uranus and Neptune as 
mostly made of ices is supporte d by the observed atm ospheric 
strong enrichment in methane dHubbard et al.l [19951) . but ice- 
free interior models are also allowed by the gravity field data 
dPodolak et al.l 120001) . We present results for the information 
content of k.2 based on our model series A-I in § 13.21 

3. Results 

3. 1 . Water in the interior of GJ 436b 

We study the effect of surface temperature, atmospheric tem- 
perature profile, and mean density on the bulk composition in 
Fig.Q] Each panel portrays the mass distribution of the three bulk 
components H/He, water and rocks. For the 3L models, panels 
show the composition in dependence on the transition pressure 
Pi _2 between the H/He envelope and the water envelope. Layers 
are separated. The water layer is subdivided in up to five re- 
gions indicating, in accordance with the water phase diagram, 
the fluid molecular and fluid dissociated phase, the dissociated 
phase (only series C, F, G), the plasma phase, the superionic 
phase, and a thin region around the first-order phase transition 
between plasma and superionic water. Due to the finite resolu- 
tion of the water phase diagram, these water layer regions are 
accurate within 0.5M e . For the 2L models (series C, F, I) panels 
show the composition in dependence on envelope metallicity Z\ . 

For all series, models without water are possible as well 
as models without rocky core, constituting the limiting cases 
of each series. Isothermal surface layers extendig down to 1 
kbar reduce the core temperature by a factor of 2 to 3 com- 
pared to the fully adiabatic case (compare e.g. series A, where 
T core = 9900 - 15000 K and B, where T core = 2700 - 6100 K). 
The same effect is observed from the assumed uncertainty in 
T\ (compare e.g. series A and G). Pressures at the core-mantle 
boundary range from 1 Mbar for large cores up to 15 Mbar for 
small cores, and the water density increases up to 6.5 g cirT 3 . 

For the fiducial M-R values, the coldest models obtained 
have deep internal temperatures as low as » 3000 K at 2-4 Mbar. 
According to the FT-DFT-MD-based H2O phase diagram, wa- 
ter then adopts the superionic phase. This phase differs from 
water-ice phases VII and X occurring at temperatures below 
2500 K at these pressures, in that protons are not confined to 
lattice positions but diffuse through the bc c oxygen lattice mak - 
ing the medium electrically conductive (iRedmer et al.l 120101) . 
Despite the low minimum surface temperature of 300 K con- 
sidered, we find no model where water is in an ice phase. The 
melting temperature is even lower for alloys of water ice mixed 



with other species at any given pressure, supporting our con- 
clusion that GJ 436b does not contain ice. Indeed, models with 
fully adiabatic profiles and those with T\ » 700 K only have 
a 0-3M ffl thin mass shell where water is fluid before it transits 
smoothly to the plasma phase because of rising temperature. For 
T\ = 700 + 100 K down to 1 kbar (e.g. series E) it is sen- 
sitive to the extension of the H/He envelope, if deep interior 
pressures become high enough for superionic or plasma water 
to be the preferred phase. Eye-striking differences between the 
3L models with and without isothermal atmosphere indicate that 
the temperature at * 1 kbar is an important quantity demanding 
accurate modeling of the atmosphere in order to narrow down 
the set of solutions. Series A models have Tiooo = 1950 K, 
series D models have Tiooo = 4000 K, and series G models 
have Tiooo = 5500 K. The value of Tiooo strongly influences 
the amount of H/He required to match R p : Mjj/jj e /M p decreases 
from 1 1-25% in series B (T woo = 300 K) to 2-11% in series H 
(r 100 o = 1300 K), further to 0.7-10% in series A to 0.1-4% in 
series D, and finally to 0.01-1% in series G. 

Since the scale heig ht of the atmosphere de creases with 
mean molecular weight ( iMiller-Ricci et al.l 120091) . outer H/He 
envelopes occupy a large volume and thus require a relatively 
low mass to account for the observed radius of planets. However, 
H and He are compressible materials so that, if mixed into the 
deep envelope where pressure rises up to 15 Mbar, the scale 
height of an H/He-water layer becomes less than that of sep- 
arated layers. Therefore, 2L models generally require a higher 
H/He mass fraction at same boundary conditions. In series C, 
the H/He mass fraction is 20.6-35.7%, in series F 13.8-19.6%, 
and in series I 5.6-14.5%, thus more than in the comparable 
cases B, E, and H. 

Increasing the mean density within the observational error 
bars of M p and R p has a neglegible effect on bulk composition 
and internal P — T profile, compare series J with E. On the other 
hand, the minimum mean density p p = 0.45 g cirT 3 increases 
the H/He mass fraction from 5-18% to 25-37%, compare series 
K and E. Furthermore, less dense models have slightly lower 
temperatures and pressures at given mass shell so that, for the 
atmospheric temperature chosen, water in the deep interior will 
be superionic in all series K models. 

The bulk composition of all models presented here is 
Mu/uJMp = 10- 4 - 0.37, M H2 o/M p = - 0.9999, and M rocks /M p 
— - 0.99. While the H/He mass fraction in models with Tiooo ^ 
4000 K can be as little as 10~ 4 - 10~ 3 , the H/He layer in such 
cases extends over 0.07 to 0.17 R p (0.30 to 0.72fl e ). We con- 
clude that H/He is an ind ispensable component of this planet. 
This is unlike GJ 1214b dCharbonneau et all 1201 Oh . where the 
measured radius ca n be explained by the as sumption of water 
steam atmospheres dRogers & Seagerll2010bl) . 

Each model series can be compared t o composition predic- 
tions from the planet formation models bv lFigueira et alj d2009l) . 
who investigate the composition of a GJ 436b-sized planet with 
separated layers of rocks+iron, ices, and H/He similar to our 
3L models. Their conditions are (2) M roc k s IM p = 0.45 - 0.7, 
(if) M H/H JM p = 0.1 - 0.2, and (Hi) M mo /M p = 0.17 - 0.40. 
Considering water representative of a mixture of elements with 
a mean density distribution similar to that of water, i.e. an ap- 
propriate mixture of ices, rocks, and H/He, these conditions 
would relax to (?) M rocks /M p < 0.7, (if) M H /He < 0.2, and 
(Hi') Mu2o/M p > 0.17. We compare in the following with the 
stricter conditions having in mind that more models might be 
possible if these conditions are relaxed. Apparently, no model 
of series A, D, and G satisfies (i - Hi), since high temperatures 
in the adiabatic outer envelopes act to decrease the mass den- 
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Fig. 1. Mass distribution of the three components H/He (yellow), water, and rocks (gray) for the 1 1 series A-K described in TableQ] 
Each panel shows various models running vertically, which differ in transition pressure Pj_2 in the case of 3L models, and in 
envelope water mass fraction Z\ for 2L models (series C, F, and I). Thick black lines show layer boundaries, and numbers close to 
filled black circles denote temperature in K and pressure in Mbar. Water is subdivided into five colors coding fluid molecular and 
fluid dissociated water (green), plasma (red), superionic water (blue), an undissolved region in the water phase diagram close to the 
superionic phase (violet), and fluid dissociated water (orange) in series C, F, and G. In series C, F and I, the H/He mass fraction can 
be obtained via (1 -Zi)x(l - M core jM p ). In series G, the center of the planet is at the top of the figure while for all other series, the 
center is at the bottom. 



sity at given pressure level, and the pressure distribution with 
mass is predominantly a function of M p (R p ) alone. As a result, 
these models have thin H/He envelopes and large cores to ac- 
count for M p (R p ). In the coldest model series B, we see agree- 
ment with (i - Hi) if Pi_2 = 115- 160 GPa. Such models have 
completely superionic water layers. In series E, we see agree- 
ment if Pi_2 = 75 - 120 GPa. Water forms an ionic fluid in a 
thin shell before it transits to the superioni c phase. This is sim - 
ilar to the interior of Uranus and Neptune (iRedmer et al.ll2010l) . 
whereas the rocky core mass fraction of 0.45-0.6 in these mod- 
els is much greater than in Uranus (< 0.18) or Neptune (< 0.24). 

In series H, we see agreement only if P1-2 ~ 80 GPa. 
However, the water shell is in the plasma phase, so that strict 
separation from the H/He envelope might not be realized. We 
conclude that models of series H are unlikely. The same argu- 
ment holds for series A, D, and G. With series C and F, strict 
and relaxed conditions give the same possible range of mod- 
els. We see agreement, if Z\ « 0.6 (F: 0.4-0.5), while the 
H/He content becomes too large for higher envelope metallici- 
ties, and the total water mass fraction too small for lower metal- 
licities. Models of series C are cold enough for water to be su- 



perionic. Whether a lattice structure of the oxygen subsystem 
or homogeneity throughout such an envelope can be maintained 
remains to be investigated. With series I, we see agreement if 
Zi = 0.55 - 0.7. These models have T core « 13,300 K and 
Pcore = 3-5.3 Mbar, where water is in a plasma state, hence 
likely to be mixed with the H/He component. The same holds for 
the allowed series F models, where T core = 8100 - 8500 K. Thus 
these models can be considered self-consistently. With series J, 
we see agreement if Pi_2 = 65 - 1 10 Mbar, where water close to 
the core is superionic justifying the assumption of a water shell. 
The fluid fraction of the water layer might be mixed into the 
H/He envelope. The same holds for series E(p = 1.69 gem" 3 ), 
which differs from series J (p = 2.22 g cuT 3 ) only in mean mass 
density p. The similarity of the range of series E and J models 
implies that the upper limit of this observational quantity is suf- 
ficiently well known. 

For seri es K models, the H/H e content is clearly more than 
predicted bv lFigueira et al.l d2009l) . since they investigate the for- 
mation of a planet with R p < 4.4R B , whereas in series K, R p = 
5.35/? e (p = 0.45 g cm' 3 ), requiring a larger fr action of light 
elements. However. iLecavelier des Etangsl d2007h predicts com- 
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plete evaporation of GJ 436b wi thin 5 Ga, if p„ < .74 g cm -3 . 
While contrary positions exist (lErkaev et al.ll2007l) . theoretical 
mass loss rates offer a way to tighten the uncertainty in age and 
maximal radius of GJ 436b. 

According to the models by ISpiegel et al.l d2010h . the at- 
mosphere is isothermal at T — 1300 K between about 1 and 
100-1000 bar, corresponding to our series H and I. From the 
discussion above, inconsistency with the water phase diagram of 
series H indicates that water is at least partially mixed within the 
H/He envelope. The extreme case of complete mixing between 
H/He and H2O is realized in series I. Models of that series with 
Zi = 0.55 - 0.70 give the best overall agreement with avail- 
able constraints. Those models have an H/He mass fraction of 
0.13-0.145, while core mass and water mass fraction cover the 
range proposed bv lFigueira et al.l (|2009l) . Alternative constraints 
may become available from improved formation models and at- 
mosphere models in the future and be applied to our broad set 
of series so that improved estimates of the bulk composition can 
then be derived from Fig. Q] 

Our models comprise those by iRogers & Seagerl d2010al) . 
who find an H/He mass fraction of 2-14.5%. Temperature- 
effects of our water equation of state tend to enhance the volume 
of a mass shell at given pressure level, requiring less H/He to 
be added to reproduce the large radius. At typical internal con- 
ditions P ~ 2 Mbar and T ~ 6000 K in GJ 436b, the tempera- 
ture eff ects enhance the vo lume by as 20% compared to water at 
300 K dFrench et alJl2009h . For this reason, we are able to find 
models with very low H/He mass fractions (« 10~ 3 M ; ,). 



3.2. Love number, core mass, and metallicity of GJ 436b 

For various models we calculated the Love number k 2 . This 
parameter is known to measure the central condensation or, 
equi valently, homogeneity of ma tter in a nearly spherical ob- 
ject ([Zharkov & Trubitsvnl 1 1978b with an upper limit of 3/2 
for a sphere of constant density and a lower limit of zero for 
Roche-type models . Jupit er's theoretical value is k 2 j ~ 0.49 
dRagozzine & WolfJ l2009h indicating a close to « = 1 poly- 
tropic density profile, which is confirmed by our calculations. 
For Neptune, we find k 2y ^ = 0.16 indicating stronger central 
condensation, for Saturn we find intermediate values around 
k 2 $ = 0.31, and for an artificial M p = 20M ffl water planet we 
find &2,w = 0.89, in agreement with water being less compress- 
ible than the H/He dominated envelopes of the solar giant plan- 
ets. Finally for GJ 436b, we find k 2 ranging from 0.02 to 0.82. 
These results are displayed in Fig. [2] together with the uncer- 
tainty in core mass of these planets. 

Models of GJ 436b with £2 = 0.02 - 0.2 are accompanied by 
any core mass between and « 0.9M p : a measurement of £2 in 
this regime would not help to constrain the interior any further. 
All our calculated 3L models with metal-free or low-metallicity 
H/He envelopes fall into this degenerate regime. Beyond the re- 
gion of strong degeneracy at k 2 > 0.2, the uncertainty in M core 
for an observationally given £2 shrinks with ki, so that upper lim- 
its for the uncertainty in core mass can be derived, if k 2 and the 
atmospheric temperature profile are known. Those models are 
necessarily two-layer models with one homogeneous envelope. 
Assuming for instance an atmospheric temperature of 1300 K, 
then any redistribution of metals from the outer part of the en- 
velope to the inner part of the envelope would (if the central 
condensation is to be kept at constant k 2 value) require a com- 
pensating transport of matter from the core upward, making the 
core mass decrease and drop below the upper limit for that at- 
mospheric temperature at any given hi. This upper limit in core 
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Fig. 2. Love number k 2 and core mass of two-layer models 
(squares connected by dashed lines), three-layer models with 
metal-free envelopes (triangles), and three-layer models with 
low-metallicity envelopes (diamonds) of GJ 436b. Color-coded 
circles are for interior models of Jupiter (red), Saturn (yellow), 
and Neptune (blue), and a 20M ffl water planet (black). Adding 
the mass of the water layer to the mass of the rocky core shifts 3L 
models solutions from filled to open triangles, and the Neptune 
model to the position as indicated by an arrow. For each of the 
dashed lines (which are for different isothermal atmospheres of 
300 K (blue), 700 K (cyan), and 1300 K (orange)), the squares 
indicate solutions for Z\ = 0, 0.5, and the highest possible value 
when M core = 0. On the left hand side of the gray dashed line 
at k 2 = 0.2, the models are highly degenerate, and on the right 
hand side, no metal-free envelope models were found. 



mass increases with decreasing atmospheric temperature, since 
dense envelopes act to reduce the density gap between envelope 
and core, which would reduce the central cendensation (£2), if it 
was not compensated for by a more massive core. If we adopt 
300 K as coldest possible temperature of an isothermal atmo- 
sphere, then the area above the blue dashed curve in Figf2]is not 
accessible by present GJ 436b. The high-core mass end of that 
curve corresponds to a zero-metallicity envelope. At the high- 
core mass end, however, the solutions lie in the area of large 
degeneracy, where besides the mass the mean density of the core 
strongly affects hi, and the behavior of solutions is no longer a 
simple function of envelope metallicity and temperature. 

At Saturn's theoretical k 2 s = 0.31, the uncertainty in core 
mass is « 0.4M p , it is » 0.2M P at Jupiter's k 2J = 0.49, and 
it is < Q.\5M p for k 2 > 0.6. Such high k 2 values, if measured, 
imply high outer envelope metallicities for this planet. For 2L 
models, we find 0.55 < Z\ < 0.78 in 300 K cold atmospheres, 
and 0.7 < Zi < 0.884 in 1300 K hot atmospheres, where k 2 
is a sensitive function of Z\ . An observational k 2 value between 
0.3 and 0.8 would help constrain envelope metallicity and core 
mass. 

Furthermore, a superionic water layer can be regarded in the 
traditional manner as a former ice layer of the initial core that 
warmed up during formation through mass accretion and con- 
traction. Adding the mass of the water layer to the underlying 
rocky core results in significantly higher core masses. The case 
of Neptune illustrates the general difficulty of applying the label 
core mass to metal-rich planets with a layered structure. 
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3.3. M-R relations for water planets 

Figure [3] shows M-R relations of water-planets with a small 
(0.25M e ) rocky core, using H 2 0-REOS. We present the M-R 
relation for profiles with an isothermal atmosphere down to 1 
kbar of 1000 K, 2000 K, for fully adiabatic profiles starting 
with ITi — 1000 K, and , for comparison, for the fit-formula from 
iFortnev et al.1 d2007bl) for pure adiabatic water planets based on a 
finite temperature correction to a zero-K isotherm of water. Our 
M-R relations for T\ = 1000 K can be fitted with the formula 



R(M) = fl () + a\ * log M + a 2 * log M" + a 3 * log M 



(1) 



with fl =1.6586, «i=0.9950, a 2 =0.1549, a 3 =0 for the solid line 
and a =2.8210, fli =-0.2928, a 2 =0.9037, a 3 =-0.1760 for the 
dashed-dot line in Fig. [3] This uncertainty in the atmospheric 
temperature profile (isothermal or adiabatic) obviously influ- 
ences the radius by ~ 0.5R®. 
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Fig. 3. Mass-radius relationships of water planets in comparison 
with Uranus (U), Neptune (N), the three known Neptune-mass 
planets, and the two known super-Earth mass planets GJ 1214b 
and CoRoT-7b. Solid line: Ti^lOOO K, isothermal atmosphere, 
Mcore =0.25M e ; dashed line: same as solid line but ri=2000 K; 
dot-dashed line: same as solid line but adiabati c atmosphere; 
thin sol id line: fit-formula for water planets from IFortnev et all 
(I2007bl) . The thin dashed-dot line is the theoretical M-R relation 
for Mg-silicate planets from I Valencia et ail (l2010h . 



H/He mass fraction < 0.01% are predicted by evolution calcula- 
tions that i nclude the energy-lim ited mass loss caused by stellar 
irradiation ( IValencia et alJl201QI) . 

4. Discussion 

4.1. A secondary planet 

To explain the high eccentricity e = 0.15 of G J 436b, the 
prese nce of a secondary planet has been sugg ested ( Ribas et all 
12.0081) but not yet confirmed by measurements dBean et alj2008l) . 
While observations exclude a resonant perturber in 2:1 reso- 
nance down to lunar mass, the presence of a secular perturber at 
distant orbit in a psidal alingment with GJ 436b is still an appeal- 
ing possibili ty (iBatvgin et al.l 2009al) . If detected, the scheme 
presented in IBatvgin et al.l d2009bl) allows determin ation of £ 2 , 
if the planets are on coplanar orbits dMardlind2010l) . We predict 
0.02 <k 2 < 0.82. 

If future measurements deny there is a secular perturber, a 
(minimum) Love number can be estimated by requiring suf- 
ficiently weak tidal interactions between planet and star that 
would allow for the planet's presence at observed orbital dis- 
tance and eccentricity. Using the common lo w-eccentricity 
approx imation of the tidal evolution equations, iJackson et al.l 
d2008l) find that an effective constant tidal Q value of 10 6 5 
results in an initial «o - <?o distribution that best matches the 
observed distribution of tidally unevolved planets. With Q' = 
3Q p /2k2 (i n the approximation of in co mpressible h omogeneous 
matter; see lGoldreich & Sotedl9^ and lLovell91 ll §123) where 
Q~ p is the tidal dissipation function, and with Q p = = 
(0.9 - 33) x 10 4 dZhang & Hamiltonll2008t iBanfield & Murray! 
1 19921) as proposed for Neptune, this implies ki = 0.0043-0.156. 
We find £ 2 > 0.02. This raises the possibility that the Q value of 
GJ 436b is close to Neptune's upper limit, making the eccentric- 
ity a consequence of weak t idal interactio n with the star. 

However, according to IJackson et al.l d2008l) the initial ec- 
centricity of GJ 436b would have been eo ~ 0.6, which has 
been demonstrated to lie beyond the range of applicability of the 
low-eccentricity approximation, as the theoretical tidal evolution 
timescale can decrease by several order s of magnitude when the 
exact orbital evolution theory is used dLeconte et al.l l2.Q10b . A 
re-calculation of GJ 436b's tidal evolution would therefore be 
helpful in order to obtain present-time (Q', k 2 ) pairs that repro- 
duce the observed orbital parameters, thereby clarifying the need 
for an external perturbing body. 



In the 5-50 M e planet mass region, previous and current M- 
R relations are similar with a flatter radius dependence in the 
latter case indicating a smaller polytropic exponent y. In par- 
ticular, the interior profiles of our water planets can be approx- 
imated well by y — 2.8 for 2 < p < 4 g cirT 3 and y - 2.4 
for p > 4 g cm 4 . The resulting theoretical radii of water plan- 
ets are below the observed radius of the thr ee currently know n 
Neptune-mass planets GJ 436b, HA T-P-1 lb dBakos et al.ll2010l) . 
and Kepler-4b dBorucki et al.ll2010h . indicating the presence of 
low-Z material in these planets. In case of warm atmospheric 
temperatures of 1000 K or more, an outer H/He envelope that 
contributes in radius ~ 15% may contribute <sc 1% in mass. The 
only known planet whose observational parameters M p , R p , and 
T eq can be explained by the assumpti on of a water-dominated 
water +rock composion without H/He dRogers & Seagerll2010bl) 
is GJ 1214b, whereas for the super-Earth CoRoT-7b dLeger et al.l 
120091) . M p = 4.8 + 0.8, a water mass fraction < 10% and an 



4.2. Constraints from formation and atmosphere profiles 

Among the broad sets of models presented here, those that 
are consistent with available constraints from formation theory, 
model atmospheres, and solubility of fluid water in H/He have 
a warm outer envelope with water mass fraction 55-70%, and a 
rocky core. At such large envelope heavy element abundances, 
ki begins to be sensitive to core mass. On the other hand, mea- 
suring fc 2 > 0.3 would either imply a violation of the condition 
from formation theory M rnc klM p > 0.45 or indicate the presence 
of rocks in the envelopes (partially replacing water). In the outer 
part of the envelope, water is in a supercritical fluid state and then 
transits to the plasma phase because of rising temperatures. This 
supports the assumption that water is mixed into the H/He enve- 
lope in the two layer models instead of forming a separate layer 
as in the three layer models. Depending on the process of forma- 
tion, however, a separation into few distinct layers with different 
water fractions can be possible. Colder models where water is 
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in a superionic state at 3000-6000 K are not consi s tent w ith the 
as 1300 K hot model atmosp heres bv | Spiegel et all (1201 Ol) or the 
recent model atmospheres bv lLewis et al.l(l2010l) which predict a 
radiative layer between 1 and 100 bar at 1 100 to 1200 K depend- 
ing on metallicity. High-pressure ice phases occur at even lower 
temperatures, hence in none of our models. Considering water 
as a proxy for metals, higher H/He mass fractions are possible 
than the 0.13-0.145 of the preferred models. The broad set of 
models presented here demonstrates the importance of accurate 
model atmospheres. Neglecting the available constraints from 
formation and atmosphere models, the set of solutions satisfying 
only M p , R p and consistency with T pe ff includes models with 
M H /hJM p = 10- 3 and M H2 o/M p = 0.999. 

4.3. Constraints from atmosphere observations 

For GJ 436b, an observational value > 0.5 implies a rocky 
core mass below 0.2M p and an envelope metallicity over 50%. 
An observational value > 0.2 would imply M core < Q.55M p , 
and Zi above ~ 40% if the envelope is homogeneous. While ho- 
mogeneous envelope models are supported by the hot internal 
temperatures and miscibility of water and hydrogen, we cannot 
be sure about the absence of envelope layer boundaries, so can 
only obtain upper limits for the core mass for any given Ic2 value. 
In order to break this degeneracy further, an estimate of the outer 
envelope metallicity is required. 

For solar system giant planets, atmospheric heavy element 
abundances are only known within a factor of a few for some 
species that do not condense at optically thick deeper levels. 
The abundances are believed to be representative of the outer 
envelope because of vertical mixing still at high altitudes be- 
low 1 bar. Irradiated extrasolar planets, such as GJ 436b how- 
ever, become optically thick at pressures several orders of mag- 
nitude lower than the onset of the convective envelope. Even in 
the fortunate case of measured atmospheric abundances, rep- 
resentative metallicity determinations might require including 
non-equilibrium chemistry, cloud formation, vertical mixing be- 
yond mixing-length theory, and other effects. On the other hand, 
current theoretical and observational progress is encouraging. 
Transmission and emission spectra predicted by state-of-the-art 
model atmospheres exhibit various observationa l signatures in 
dependence on metallicity. IChabrier et al.l ( |2007|) found a 30% 
decrease of a hot Jupiter's emission at 4.5pm due to encanced ab- 
sorption by CO when the overall metallicity was enhanced from 
lx to 5x solar. For super-Earth atmospheres, iMiller-Ricci et al.l 
(120091) found large hydrogen fractions (low mean molecular 
weight) to map on a detectable variability of the transit radius 
with wavelength and a rapid decrease in va riability with mean 
molecular weight. For GJ 436b in particular. iLewis et al.l d2010h 
developed a three-dimensional atmospheric circulation model 
and predict metal-enhanced atmospheres to reveal themselves 
in an orbital phase dependency of spectra potentially detectable 
during future missions. In ad dition, secondary eclipse measure- 
ments (IStevenson et al.ll2010l) already provi de hints of a strong ly 
(30-50x solar) metal-enhanced atmosphere (ILewis et a l. 2010). 

4.4. Comparison with Uranus and Neptune 

Interior models of Uranus (M p = 14.5M e ) predict that at least 
21% by mass is H/He, and at least 1 5% in Neptune (M p = 
17.1M e ) (iFortnev & Nettelmannll2010l) . This is consistent with 
the preferred model values for GJ 436b (13-14.5%) considering 
that some fraction of water as assumed in the model calculations 



could indeed be rocks and H/He in the envelope of either of the 
planets. On the other hand, since only lower limits of the H/He to 
ice and rock ratios in these planets can be derived, their compo- 
sition can also be very different from each other. Irradiation by 
the close-by parent star may rise the core temperatures of pre- 
ferred GJ 436b models significantly above those of Uranus and 
Neptune, where adiabatic models predict T core = 6000 - 7000 K, 
hence water in the superionic phase. Methane has been sug- 
gested as separat ing into metallic hydrogen and diamond (e.g. 
iHirai et aT]l2009l) at those deep interior conditions, so that the 
interior of those planets might be organized very different by 
that from GJ 436b, even in the case of similar bulk composition. 
Interior models of Uranus and Neptune require supersolar heavy 
element abundance in the outer envelope and some admixture of 
H/He in the deep interior, which can be represented by at least 
two envelopes with different compositions. Such a structure is 
bracketed by our 2L and 3L models. The observed low heat flux 
of Uranus and, to a lesser extend, of Neptune gives further evi- 
dence of layer boudaries, or more general, inhomogeneities in 
these planets dHubbard et al.l [l995l) . Unfortunately, the age of 
GJ 436 is poorly known, but investigations of the cooling be- 
havior of various interior models can help further rule out or 
strenghten some of the models presented here. 

4.5. Conclusions 

Applying H-, He-, and H2O-REOS for hydrogen, helium, and 
water to generate three-layer models (H/He envelope, water en- 
velope, rocky core) and two-layer models (H/He/HiO envelope, 
rocky core) of GJ 436b, we find interior models ranging from 
large rocky cores with significant H/He layers to almost pure 
water worlds with thin, but extended H/He atmospheres. If water 
occurs and the 1 kbar temperature is not much higher than 700 K, 
this exoplanet constitutes the third natural potential object for 
realyzing superionic water, besides the deep interior of Uranus 
and Neptune. In the H/He envelope, the temperature rises rapidly 
such that we do not find high-pressure ice phases in the inner en- 
velope. Indeed, if we assume temperature profiles in good agree- 
ment with model atmospheres, water in the interior of GJ 436b 
transits from fluid to plasma, favoring an outer envolope with 
strong water enrichment instead of a separation into two distinct 
envelopes. 

The Love number Ic2 of different possible compositions with 
metal-free outer H/He envelopes varies from 0.02 to 0.2. In this 
£2 regime, the core mass depends non-uniquely on £2 and can- 
not be derived from a measurement of £2- An observationally 
derived value k 2 > 0.2 (0.3; 0.5) would imply M core < Q.55M p 
(M core < 0.4M p ; M core < 0.2M p ) and a high outer envelope 
metallicity. 

The radius of GJ 436b is larger than predicted for warm- 
water planets. Even in case of the warmest models found, a small 
H/He fraction of 10~ 3 M ; , is necessary to account for the radius. 
Calculating the cooling behavior on the basis of a proper treat- 
ment of the atmosphere might help to better constrain the interior 
of this planet. 

With this work we want to motivate both accurate transit- 
timing observations of this planet that aim to clarify the exis- 
tence of a secondary planet and that may allow derivation of 
kj_ and transmission spectra measurements aimed at determining 
the atmospheric metallicity. 
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